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PREFACE 


Welcome to Control System Design with LabVIEW!M 


As a companion to the textbook, Modern Control Systems by Richard C. Dorf and Robert H. Bishop, this supplement provides a 
set of comprehensive tutorials and exercises utilizing the LabVIEW Control Design and Simulation Module. NI LabVIEW is a 
graphical and textual language for prototyping real control and signal processing systems enabling a seamless flow from simulation 
to deployment on real-time control hardware. This text guides students through the process of modeling a system, analyzing the 
model to build a controller, and validating the robustness of the control system thorough simulation. Its contents and organization 
mirror the corresponding sections in MCS making it an ideal teaching companion, especially when LabVIEW is being utilized in a 
corresponding controls laboratory. 

It is assumed readers will have access to LabVIEW 2010 or later, the Control Design and Simulation Module, and Mathscript 
RT. With these tools the reader can easily build, simulate, and analyze the examples and problems included in the text. Solution VIs 
for many of the problems are included for use as a programming reference and can also be used as a starting point for solving more 
advanced design problems. All of the LabVIEW examples were developed and tested on a PC compatible with LabVIEW Express 
2010. The available VIs can be downloaded from http://www.ntspress.com/publications/modern-control-systems-with-labview. 

For students unfamiliar with LabVIEW, it will be very helpful to have access to the Learning with LabVIEW textbook by Robert 
H. Bishop, available from Prentice Hall. For readers new to LabVIEW control, a wealth of documentation exists. Please see 
www.ni.com/academic/controls.htm. 

We wish to express appreciation to Andy Chang at National Instruments and to Jorge Alvarez at The University of Texas at Austin 
for their work in creating the LabVIEW VIs and the screen captures used in the making of this companion text. Special thanks also 
to the folks at National Instruments, especially Erik Luther and Dr. Jeannie Falcon, for their continued support. 


ROBERT H. BISHOP 


Mathematical Models of Systems 


Application of the many classical and modern control system design and analysis tools is based on mathematical models. LabVIEW 
can be used with systems given in the form of transfer function descriptions. We begin this chapter by showing how to use LabVIEW 
to assist in the analysis of a typical spring-mass-damper mathematical model of a mechanical system. Using a LabVIEW virtual 
instrument (or VI) we can develop an interactive analysis capability to analyze the effects of natural frequency and damping on the 
unforced response of the mass displacement. We also discuss transfer functions and block diagrams. In particular, we are interested 
in how LabVIEW can assist us in manipulating polynomials, computing poles and zeros of transfer functions, computing closed-loop 
transfer functions, computing block diagram reductions, and computing the response of a system to a unit step input. The chapter 
concludes with a design example for an electric traction motor control design. 


1.1 Simulating Spring-Mass-Damper Systems 


A spring-mass-damper mechanical system is shown in Fig. 1.1. The motion of the mass, denoted by y(t), is described by the 
differential equation 


MY) + by(t) +ky@) = r(t). 


The unforced dynamic response of the spring-mass-damper mechanical system is 


yO) _ . 
y(t) = ———— e **" sin (onv — ¿2? t+ e) 
y/1-z? 


where 0 = cos! ¢, c»? = k/M and 2¢a, = b/M. The initial displacement is y(0) and y(0) = 0. The transient system response is 
underdamped when ¢ < 1, overdamped when ¢ > 1, and critically damped when ¢ = 1. 


Example 1.1  Spring-Mass-Damper Simulation 


We can use LabVIEW to visualize the unforced time response of the mass displacement following an initial displacement of y (0). 
Consider the underdamped case, where 


1 
0)—0.15m, o, = V2 rad/sec, — ——. (k/M =2, b/M — V). 
y(0) m, @ rad/sec, ¢ Wp (k/ / ) 


Wall 
friction, b 





r(t) 


Force 


Figure 1.1: A mass-spring-damper system. 
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Figure 1.2: VI to analyze the spring-mass-damper. 


The LabVIEW commands to generate the plot of the unforced response are shown in Fig. 1.2. 

In the LabVIEW setup, the variables y (0), œn, £, and t are input to the user interface part of the VI. Then the Unforced.vi is 
executed to generate the desired plots. This creates an interactive analysis capability to analyze the effects of natural frequency 
and damping on the unforced response of the mass displacement. One can investigate the effects of the natural frequency and the 
damping on the time response by simply entering new values of w, and ¢ and rerun the Unforced.vi. © 


For the spring-mass-damper problem, the unforced solution to the differential equation was readily available. In general, when 
simulating closed-loop feedback control systems subject to a variety of inputs and initial conditions, it is difficult to obtain the solution 
analytically. In these cases we can use LabVIEW to compute the solutions numerically and to display the solution graphically. 


1.2 Analyzing Systems Using LabVIEW 


LabVIEW can be used to analyze systems described by transfer functions. Since the transfer function is a ratio of polynomials, we 
begin by investigating how LabVIEW handles polynomials, remembering that working with transfer functions means that both a 
numerator polynomial and a denominator polynomial must be specified. In LabVIEW, polynomials are represented by row vectors 
containing the polynomial coefficients. For example, the polynomial 


P(s) =s? +357 +2543 


is entered as [3 2 3 1] as shown in Fig. 1.3. 
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Figure 1.3: Entering the polynomial P (s) = s? + 3s? + 2s +3 and calculating the roots of P(s) = 0. 


If p is a row vector containing the coefficients of P (s) in ascending order, then Polynomial Roots (p) is a row vector containing 
the roots of the polynomial. Conversely, if r is a row vector containing the roots of the polynomial, then Create Polynomial 
from Roots(r) is a row vector with the polynomial coefficients in ascending order. We can compute the roots of the polynomial 
P(s) = s? + 3s? + 2s + 3 with the Polynomial Roots function. In Fig. 1.3, we show how to reassemble the polynomial with the 
Create Polynomial from Roots function. 

Multiplication of polynomials is accomplished with the Multiply Polynomials.vi function. Suppose we want to expand the 
polynomial N (5), where 


N(s) = Qs? +25 + 1)(s +4). 


The associated LabVIEW commands using the Multiply Polynomials.vi function are shown in Fig. 1.4. Thus, the expanded 
polynomial is 
N (s) = 3s? + 14s? + 9s + 4. 


The function Polynomial Evaluation is used to evaluate the value of a polynomial at the given value of the variable. The polynomial 
N (s) has the value N(—5) = —66, as shown in Fig. 1.4. 

The LabVIEW Control Design & Simulation Toolbox treats linear, time-invariant system models as objects, allowing you to 
manipulate the system models as single entities. In the case of transfer functions, you create the system models using the CD 
Construct Transfer Function Model; for state variable models you employ the CD Construct State Space Function Model. 
The use of CD Construct Transfer Function Model is illustrated in Fig. 1.5. 

Based on the LabVIEW object-oriented programming capabilities, the system model objects possess object properties that can 
be modified; likewise functions that operate on system model objects are called methods. For example, if you have the two system 


models 


10 1 
Gi(s) = 312545 and G2(s) = "EE 


you can add them using the Add Rational Polynomials.vi function to obtain 


s? + 125 +15 


G(s) =G G SS 
(8) = Gils) + G2) = Sa r3 
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Figure 1.4: Using Multiply Polynomials and Polynomial Evaluation to multiply and evaluate the polynomials N (s) = (3s? 4- 2s + 
1)(s 4- 4). 
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Figure 1.5: Using the CD Construct Transfer Function Model, the Add Polynomials, and the Draw Transfer Function Equations 
functions. 
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Figure 1.6: Using the CD Pole-Zero Map function to compute the poles and zeros of a linear system. 


The corresponding LabVIEW commands are shown in Fig. 1.5. Computing the poles and zeros associated with a transfer function 
is accomplished by operating on the system model object with the CD Pole-Zero Map function, as illustrated in Fig. 1.6. In the 
next example, we obtain a plot of the pole and zero locations in the complex plane. On the pole and zero map, zeros are denoted by 
an “o" and poles are denoted by a *x". 


Example 1.2 Transfer Functions 


Consider the transfer functions 


Gis) = 6s? +1 
$3 4+ 35243541 
and 
H(s) = (s + Ds + 2) 


(s + 2i)(s — 2i)(s 4-3) 


Utilizing LabVIEW, we can compute the poles and zeros of G(s), the characteristic equation of H(s), and divide G(s) by H (s). 
We can also obtain a plot of the pole-zero map of G(s)/H(s) in the complex plane. The pole-zero map of the transfer function 
G(s)/H(s) and the associated LabVIEW commands are shown in Fig. 1.7. 

The pole-zero map shows clearly the five zero locations, but it appears that there are only two poles. This cannot be the case, since 
we know that the number of poles must be greater than or equal to the number of zeros. Using the CD Pole-Zero Map function, we 
ascertain that there are in fact four poles at s = —1. Hence, multiple poles or multiple zeros at the same location cannot be discerned 
on the pole-zero map. O 


1.3 Block Diagram Models 


Suppose we have developed mathematical models in the form of transfer functions for the plant, represented by G(s), and the 
controller, represented by G (s), and possibly other system components such as sensors and actuators. Our objective is to interconnect 
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Figure 1.7: The CD Pole-Zero Map function. 


these components to form a control system. We can use LabVIEW functions to carry out the block diagram transformations. A 
simple open-loop control system can be obtained by interconnecting a plant and a controller in series as illustrated in Fig. 1.8. We 
can use LabVIEW to compute the transfer function from R(s) to Y (s), as will be illustrated in Example 1.3. 







Controller TO Process 
G(s) G(s) 


Figure 1.8: Open-loop control system (without feedback). 
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Figure 1.9: Series connection block diagram. 


Example 1.3 Series Connection 


Let the process represented by the transfer function G(s) be 





1 
G(s) = ——, 
(5) = Soo? 
and let the controller represented by the transfer function G, (s) be 
s+l 
G = . 
ce(s) s42 


Suppose we want to cascade two transfer functions G(s) and G,(s), as illustrated in Fig. 1.9. The transfer function G.(s)G(s) is 
computed using the CD Series function as shown in Fig. 1.10. 
The resulting transfer function, Ge(s)G (s), is 


s+1 
500s3 + 10005? ` 


Gc(s)G(s) = 





Numerator of G(s) 










Numerator of Gecis] 





Figure 1.10: The CD Series vi. 
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System 1 
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Figure 1.11: Parallel connection block diagram. 
© 


Block diagrams often have transfer functions in parallel as shown in Fig. 1.11. In such cases, the function CD Parallel can be 
useful. 


Example 1.4 Parallel Connection 


Suppose we desire to connect two systems G(s) and G2(s) in parallel as shown in Fig. 1.11, where 


Gi(s) = and G2(s) = —— 


1 
500s2 
The CD Parallel function is illustrated in Fig. 1.12. The resulting transfer function G1 (s)G2(s) is 
5005s? + 500s? +s +2 


016)62) = — $5553 F 10002 





T3 Parallel Connection.vi Block Diagram 
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Figure 1.12: The CD Parallel vi. 
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Figure 1.13: A basic control system with unity feedback. 


We can introduce feedback into the control system by closing the loop with unity feedback, as shown in Fig. 1.13. The signal 
E,(s) is an error signal; the signal R(s) is a reference input. In this control system, the controller is in the forward path and the 
closed-loop transfer function is 

G.(s)G 
ro- EE) _ 

1 + Gc(s)G(s) 

We can utilize the CD Feedback function to aid in the block diagram reduction process to compute closed-loop transfer functions 
for single- and multiple-loop control systems. When the closed-loop control system has unity feedback, we can use the CD Unit 
Feedback function to compute the closed-loop transfer function. 


Example 1.5 Unity Feedback 


In this example, we consider the CD Unit Feedback function. Let the process, G(s), and the controller, G, (s), be as in Fig. 1.14. 
We compute Ge(s)G (s) using the CD Unit Feedback function to close the loop. The command sequence is shown in Fig. 1.15 and 
results in the closed-loop transfer function 


Ge(s)G(s) sl 


1+ G,(s)G(s)  500s?--1000s^ +s+ 1 


Q 


Another basic feedback control configuration is shown in Fig. 1.16. In this case, the controller is located in the feedback path. 


The closed-loop transfer function is 
G(s) 


~ 14 G(s)H(s) 
The use of the CD Feedback function for a nonunity feedback system is illustrated in Example 1.6. 


T (s) 


Example 1.6  Nonunity Feedback 


Let the process, G(s), and the controller, H (s), be as in Fig. 1.17. To compute the closed-loop transfer function with the controller 
in the feedback loop we use the CD Feedback function as shown in Fig. 1.18. The closed-loop transfer function is 


s+] 


i A S A 
(9) = $0053 4 100052 be 4 | 


© 


The LabVIEW functions CD Series, CD Parallel, and CD Feedback can be used as aids in block diagram manipulations for 
multiple-loop block diagrams, as will be illustrated in Example 1.7. 


Example 1.7  Multiloop Reduction 





Figure 1.14: Feedback connection block diagram. 
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Figure 1.17: Application of the CD feedback. 
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Figure 1.18: Application of the CD feedback function for nonunity feedback. 


A multiloop feedback system is shown in Fig. 1.19. The objective is to compute the closed-loop transfer function 








Y 
noe, 
R(s) 
when 
ent ee G2(s) = 
ppt TO! AAS) = s+? 
2 
s^ +1 s+] 
G = ———————, G ———-, 
a(s) s? +4s +4 m s T6 
and 
I e s Ho(s)=2, Hy(s)=1 
1 al rr. 2 n , 3 zm: d 


For this example, a five-step procedure is followed: 
Step 1: Input the system transfer functions into LabVIEW. 
Step 2: Move H, behind G4. 
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Y(s) 





Figure 1.19: A multiloop feedback system. 


Step 3: Eliminate the G3G4H loop. 
Step 4: Eliminate the loop containing H2. 
Step 5: Eliminate the remaining loop and calculate T (s). 


The five steps are utilized in Fig. 1.20, and the corresponding block diagram reduction is shown in Fig. 1.21. 
The result of executing the LabVIEW commands is 


0.08354 + 0.255? + 0.2552 + 0.25s + 0.167 


TG)- 3 4 3 2 
$? + 16.0835^ + 72.75s? + 137s4 + 123,667s + 59.333 


If we had performed the calcuations by hand, we would have obtained 


5? +454 + 653 + 657 +55 +2 


T(s) = 6 5 4 3 2 
1259 + 2055s? + 1066s* + 25175? + 3128s- + 2196s + 712 


We must be careful in calling this the closed-loop transfer function. The transfer function is defined to the input—output relationship 
after pole-zero cancellations. If we compute the poles and zeros of T (s), we find that the numerator and denominator polynomials 
have (s + 1) as acommon factor. This must be canceled before we claim we have the closed-loop transfer function. To assist us 
in the pole-zero cancellation, we can use the CD Minimal Realization function. The CD Minimal Realization function, shown in 
Fig. 1.22, removes common pole-zero factors of a transfer function. After the application of the CD Minimal Realization function, 
we find that the order of the denominator polynomial has been reduced from six to five, implying one pole-zero cancellation. © 


Example 1.8 Electric Traction Motor Control 


In this example we consider the electric traction motor system represented by the block diagram shown in Fig. 1.23. 

The objective is to compute the closed-loop transfer function and investigate the response of c (s) to a commanded wg(s). The 
first step, as shown in Fig. 1.24, is to compute the closed-loop transfer function w(s)/@g(s) = T (s). 

The closed-loop characteristic equation is second-order with wm, = 52 and ¢ = 0.012. Since the damping is low, we expect 
the response to be highly oscillatory. We can investigate the response w(t) to a reference input, c(t), by utilizing the CD Step 
Hesponse function. The CD Step Response function calculates the unit step response of a linear system. The step response of 
the electric traction motor is shown in Fig. 1.25. As expected, the wheel velocity response, given by y(t), is highly oscillatory. Note 
that the output is y(t) = w(t). 

© 
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Figure 1.20: Reducing the multiloop feedback system. 
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Figure 1.21: Block diagram reduction of the system in Fig. 1.20. 
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Figure 1.22: Application of the CD Minimal Realization function. 
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Figure 1.23: Block diagram of an electric traction motor system. 
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Figure 1.24: Electric traction motor block reduction. 
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Figure 1.25: The CD Step Response function. 


State Variable Models 


The time-domain method utilizes a state-space representation of the system model, given by 
x = Ax+ Bu and y = Cx + Du. (2.1) 


The vector x is the state of the system, A is the constant n x n system matrix, B is the constant n x m input matrix, C is the constant 
p x n output matrix, and D is a constant p x m matrix. The number of inputs, m, and the number of outputs, p, are taken to be one, 
since we are considering only single-input, single-output (SISO) problems. Therefore y and u are not bold (matrix) variables. The 
main elements of the state-space representation in Eq. (2.1) are the state vector x, the input u, the output y, and the constant matrices 
(A, B, C, D). 

Given a transfer function, we can obtain an equivalent state-space representation and vice versa. The function CD Convert to 
State-Space Model is used to convert a transfer-function representation to a state-space representation; the function CD Convert 
to Transfer Function Model is used to convert a state-space representation to a transfer function. An analysis of a third-order system 
is presented in Example 2.1. 


Example2.1  Third-Order System 

Consider the third-order system 

Y (s) 25? + 8s +6 
R(s) — 534+ 852 4+ 16s 6^ 





T(s) = (2.2) 


We can obtain a state-space representation using the CD Convert to Sate-Space Model function, as shown in Fig. 2.1. The 
state-space representation of Eq. (2.2) is given by Eq. (2.1) where 


—8 —4 -1.5 0.25 
A= 4 0 0 |, B= 0 ; 
0 1 0 0 
and 
C=[8 8 6], D = [0]. 
The state-space representation of the transfer function in Eq. (2.2) is depicted in block diagram form in Fig. 2.2. © 


The time response of the system in Eq. (2.1) is given by the solution to the vector differential equation 
t 
x(t) = exp(At)x(0) + | exp [A(t — t)] Bu(t)dt. (2.3) 
0 


The matrix exponential function in Eq. (2.3) is the state transition matrix, ® (t), where 
(t) = exp(Ar). 


We can use the function Matrix Exp to compute the transition matrix for a given time interval, as illustrated in Fig. 2.3. An RLC 
network is investigated in Example 2.2. 
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-8 -4-15 0.25 
4 0 0 |x(t)+] 0 |u(t) 
01 0 0 


y(t) = E 8 s| x(t) + fo u(t) 


i 2 
s +85 +165+6 


dx/dt = 





Figure 2.1: Conversion of transfer functions to a state-space representation 


Example 2.2 State Transition Matrix of a RLC Network 


The state of the system can be described in terms of the state variables (x1, x2), where x, (t) is the capacitor voltage, uv, (t), and xo(t) 
is the inductor current, iz (t). The state-space representation of the RLC circuit is 


1 1 
+=] : “ext | s and y=[0 RIx-[0]u. 
L L 


where L = inductance (H), R = resistance (Q), and C = capacitance (F). When R = 3 Q, L = 1 H, and C = 1/2 F, we have 


0 —2 2 
a=] | a B=| 4 |. C-[0 3], D=[0]. 
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Figure 2.3: Computing the state transition matrix for a given time, At = dt. 
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Figure 2.4: Computing the time response for nonzero initial conditions and zero input using CD Linear Simulation. 


The initial conditions are x;(0) = 1, x2(0) = 1.25, and the input u(t) = 0. At t = 0.2, using the state transition matrix, the state is 


X1 _ | 0.9671  —0.2968 X1 {06 
( x2 ) 7 | 0.1484 0.5219 | ( X2 ) 7 ( 0.8 ) 
The time response of the system of Eq. 2.2 can also be obtained using the CD Linear Simulation function. The CD Linear 
Simulation function can accept nonzero initial conditions, as well as an arbitrary input function. Using the CD Linear Simulation 
function, we can calculate the response for the RLC network as shown in Fig. 2.4. The state at t = 0.2 is predicted with the CD 
Linear Simulation function to be x; (0.2) = 0.6 and x2(0.2) = 0.8. We can compare the results obtained by the CD Linear 
Simulation function by multiplying the initial condition state vector by the state transition matrix to find the same result. In the 


example illustrated in Fig. 2.4, we set B = 0 to obtain the zero-input response. This allows us to keep the Square Waveform 
function input in the code for future analysis by setting B back to its non-zero values. Q 


Feedback Control System Characteristics 


The advantages of feedback will be illustrated with two examples in this chapter using LabVIEW in the control system analysis. In 
the first example, we introduce feedback control to a speed tachometer system in an effort to reject disturbances. The reduction in 
system sensitivity to plant variations, adjustment of the transient response, and reduction in steady-state error will be demonstrated 
in a second example, the English Channel boring machine. 


Example 3.1 Speed Control System 


The open-loop block diagram description of the armature-controlled dc motor with a load torque disturbance, T4 (s), is shown in 
Fig. 3.1. The values for the various parameters are given in Table 3.1. We have two inputs to our system, V,(s) and Ty(s). Relying 
on the principle of superposition, which applies to our linear system, we consider each input separately. To investigate the effects 
of disturbances on the system, we let V,(s) = 0 and consider only the disturbance 75(s). Conversely, to investigate the response 
of the system to a reference input, we let Tg (s) = O and consider only the input V,(s). The closed-loop speed tachometer control 
system block diagram is shown in Fig. 3.2. 


Disturbance 
T Pi ) 


w(s) 
Speed 





Figure 3.1: Open-loop speed control system (without external feedback). 


Table 3.1: Tachometer Control System Parameters 


Ra Km J b Kp Ka Ki 
1 10 2 05 O1 54 1 


If the system displays good disturbance rejection, then we expect the disturbance Ty (s) to have a small effect on the output œ (s). 
Consider the open-loop system in Fig. 3.3 first. We can use LabVIEW to compute the transfer function from Tg (s) to c (s) and 
evaluate the output response to a unit step disturbance (that is, 7; (s) = 1/s). The time response to a unit step disturbance is shown 
in Fig. 3.3. 

The open-loop transfer function is 

w(s) —] 
G(s) = = —————, 
Ta(s) 2s +1.5 
Since the desired value of w(t) is zero (remember that V,(s) = 0), the steady-state error is just the final value of w(t), which is 
denoted by c(t). The steady-state error, shown on the plot in Fig. 3.3, is approximately the value of the speed when t = 7 seconds. 
We can obtain an approximate value of the steady-state error by looking at the plot in Fig. 3.3. The approximate steady-state value 
of e(t) is 





wol) © o (7) = —0.66 rad/s. 


22 


CHAPTER 3 FEEDBACK CONTROL SYSTEM CHARACTERISTICS 23 


Tas) 
| Amplifier | 


T (8) ~ Tue) 1 
eos s Ho-gpérm GSH 





= EN 
K, 


Figure 3.2: Closed-loop speed tachometer control system. 


In a similar fashion, we begin the closed-loop system analysis by computing the closed-loop transfer function from 7y(s) to w(s) 
and then generating the time-response of c(t) to a unit step disturbance input. The output response and the vi are shown in Fig. 3.4. 
The closed-loop transfer function from the disturbance input is 


os) —] 
Ta(s) 2s + 541.5` 





Gl = tf(1, [J b]; 
G2 = tf(&m,1); 
G3 = tf(1 Ray 
G4 = tf(Kb,1); 
G5 = tf(-1,1); 











Figure 3.3: Analysis of the open-loop speed control system. 
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Figure 3.4: Analysis of the closed-loop speed control system. 


As before, the steady-state error is just the final value of w(t), which is denoted by c(t). The steady-state error is shown on the 
plot in Fig. 3.4. We can obtain an approximate value of the steady-state error by looking at the plot in Fig. 3.4. The approximate 
steady-state value of c(t) is 

c (00) © @o(0.02) = —0.002 rad/s. 


We generally expect that |w.(0o) /@,(0o)| < 0.02. The ratio of closed-loop to the open-loop steady-state speed output due to a unit 


step disturbance input, in this example, is 


oel) _ 9.993 





Wo (00) 


We have achieved a remarkable improvement in disturbance rejection. It is clear that the addition of the negative feedback loop 
reduced the effect of the disturbance on the output.This demonstrates the disturbance rejection property of closed-loop feedback 
systems. © 
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Figure 3.5: A block diagram model of a boring machine control system. 


Example 3.2 English Channel Boring Machines 


The block diagram description of the English Channel boring machine is shown in Fig. 3.5. 
The transfer function of the output due to the two inputs is 


Y(s) — LUE D LM + 


l 
PEEKEONENE r3 
34 1254K ols pK 40 


The effect of the control gain K on the transient response is shown in Fig. 3.6 along with the vi used to generate the plots. 
Comparing the two plots, it can be seen that decreasing K decreases the percent overshoot. Although it is not as obvious from 

the plots in Fig. 3.6, it is also true that decreasing K increases the settling time. This can be verified by taking a closer look at 

the data used to generate the plots. This example demonstrates how the transient response can be altered by feedback control gain 
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Figure 3.6: The response to a step input with K = 100 and K = 20 and LabVIEW script. 
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Figure 3.7: The response to a step disturbance with K = 100 and K = 20. 


K. Based on our analysis thus far, we would prefer to use K = 20. Other considerations must be taken into account before we 
can establish the final design. Before making the final choice of K, it is important to consider the system response to a unit step 


disturbance, as shown in Fig. 3.7. 
We see that increasing K reduces the steady-state response of y(t) to the step disturbance. The steady-state value of y(t) is 


0.05 and 0.01 for K = 20 and 100, respectively. The steady-state errors, percent overshoot, and settling times (2% criteria) are 
summarized in Table 3.2. 


Table 3.2: Response of the Boring Machine Control System for K = 20 and K = 100 


K=20 kK =100 
P.O. 4% 22% 
T; 1.0s 0.7s 
ess 5% 1% 


The steady-state values are predicted from the final-value theorem for a disturbance input as follows: 


1 1 1 
lim y(t) = lim s | ———————_ |- = —. 
t— 00 s>0 | s(s+12)+ K] 5 K 
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Figure 3.8: System sensitivity to plant variations (s = jæ). 


If our only design consideration is disturbance rejection, we would prefer to use K = 100. We have just experienced a common 
trade-off situation in control system design. In this particular example increasing K leads to better disturbance rejection, whereas 
decreasing K leads to better performance (that is, less overshoot). The final decision on how to choose K rests with the designer. 
Although LabVIEW can certainly assist in the control system design, it cannot replace the engineer's decision-making capability 
and intuition. The final step in the analysis is to look at the system sensitivity to changes in the plant. The sensitivity function is 
computed as 
s(s + 1) 
s(s--12) --K- 


We can compute the values of S(s) for different values of s and generate a plot of the system sensitivity. For low frequencies, we 
can approximate the system sensitivity by 


S(s) — 


S(s) X — 


Increasing the gain K reduces the system sensitivity. The system sensitivity plots when s = jæ are shown in Fig. 3.8 for K = 20. 
© 
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Performance of Feedback Control Systems 


In this chapter we begin our investigation of the performance of feedback control systems. As an illustrative example, we solve the 
problem of bank angle control of an aircraft using LabVIEW to build a reasonably simple yet effective model of the response to 
aileron deflections. Then using the idea of model simplification, we present a second-order approximate model. The second-order 
approximation gives us insight into the expected behavior of the system and how to obtain an initial design. Finally we use a 
LabVIEW simulation to perform analysis in an interactive mode. We begin by investigating time-domain performance specifications 
given in terms of transient response to a given input signal and the resulting steady-state tracking errors. 


4.1 Time-Domain Specifications 


Time-domain performance specifications are generally given in terms of the transient response of a system to a given input signal. 
Because the actual input signals are generally unknown, a standard test input signal is used. Consider the second-order system shown 
in Fig. 4.1. 


Plant, G(s) 


sg? + CWS 





Figure 4.1: Single-loop second-order feedback system. 


The closed-loop output is 


2 


Co 
Y (s) = —————————— R(s). 
(s) s? + 2t£cons + c? ‘s) 


We have already discussed the use of the step function to compute the step response of a system. Now we address another important 
test signal: the impulse function. The impulse response is the time derivative of the step response. The response of a second-order 
system to a step function is shown in Fig. 4.2 for various values of œn and ¢. 


Similarly, the response of a second-order system to an impulse function is shown in Fig. 4.3. In the script, we set cy, = 1, which 
is equivalent to computing the step response versus «yf. This gives us a more general plot valid for any œn > O. 


We studied the CD Linear Simulation function in Chapter 2 for use with state-variable models; and now we consider its use 
with transfer function models as shown in Example 4.1. 
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Figure 4.2: Response of a second-order system to a step function. 


Example 4.1 Mobile Robot Steering Control 


The block diagram for a steering control system for a mobile robot is shown in Fig. 4.4. Suppose the steering controller, Ge(s), is 
Kı 
Ge(s) = Kp + P3 , 


When the input is a ramp, the steady-state error is 


ess = 77-7, (4.1) 
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Figure 4.3: Response of a second-order system to an impulse function. 


where Ky = KrK. The effect of the controller constant, Kz, on the steady-state error is evident from Eq. (4.1). Whenever Ky 
is large, the steady-state error is small. We can simulate the closed-loop system response to a ramp input using the CD Linear 
Simulation function. The controller gains Kp, Kz, and the system gain K can be represented symbolically in the block diagram so 
that various values can be selected and simulated. The results are shown in Fig. 4.5 for Kp = K = 1, Kr = 2, and t = 1/10. © 
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Figure 4.4: Steering control system block diagram. 


4.2 Simplification of Linear Systems 


It may be possible to develop a lower-order approximate model that closely matches the input-output response of a high-order model. 
We can use LabVIEW to compare the approximate model to the actual model, as illustrated in Example 4.2. 
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Figure 4.5: Transient response of the mobile robot steering control system to a ramp input. 
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Example 4.2 Simplified Model 


Consider the third-order system 


6 
Hes s3 -F6s2 -- 11s +6 
A second-order approximation is 
L(s) = 1.60 | 
s? + 2.584s + 1.60 
A comparison of their respective step responses is given in Fig. 4.6. © 


Example 4.3 Aircraft Roll Control using LabVIEW 


Each time we fly on a commercial airliner, we experience first-hand the benefits of automatic control systems. These systems assist 
pilots by providing pilot relief during extended flights and by improving the handling qualities of the aircraft over a wide range of 
flight conditions. The special relationship between flight and controls began in the early work of the Wright brothers. Using wind 


Denominator of H jr) | "e 
e SISO | e -— Function His) 


Numerator of L(s) cp Construct Transfer Function Model.vi 
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Figure 4.6: Step response comparison for an approximate transfer function versus the actual transfer function. 
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Figure 4.7: (a) Control of the bank angle of an airplane using differential deflections of the ailerons. (b) Bank angle control autopilot. 


tunnels the Wright brothers applied systematic design techniques to make their dream of powered flight a reality. This systematic 
approach to design was partially responsible for their success. Another significant aspect of their approach was their emphasis on 
flight controls; the brothers insisted that their aircraft be pilot controlled. Observing birds control their rolling motion by twisting 
their wings, the Wright brothers built aircraft with mechanical mechanisms that twisted their airplane wings. Today we no longer 
use wing warping as a mechanism for performing a roll maneuver, instead we control rolling motion using ailerons, as shown in 
Fig. 4.7. The Wright brothers also used elevators (located forward) for longitudinal control (pitch motion) and rudders for lateral 
control (yaw motion). Today’s aircraft still use both elevators and rudders, although the elevators are generally located on the tail 
(rearward). 

The first controlled, powered, unassisted take-off flight occurred in 1903 with the Wright Flyer I (a.k.a. Kitty Hawk). The first 
practical airplane, the Flyer III could fly figure eights and stay aloft for half an hour. Three-axis flight control was a major (and often 
overlooked) contribution of the Wright brothers. A concise historical perspective is presented in Stevens and Lewis!. 

For our aircraft we select e; = 1.4 and k = 11.4, where the aircraft dynamics are given by G(s) = JD The associated 
time-constant of the roll subsidence is t = 1/e, = 0.7 second. These values represent a fairly fast rolling motion response typical 
of an agile aircraft. For the aileron actuator model, we typically use a simple first-order system model, 


a (s) _ P 
es) sp 








(4.2) 


where e(s) = ga (s) — o(s). In this case we select p = 10. This corresponds to a time-constant of c = 1/p = 0.1 second. This is a 
typical value consistent with a fast response. We need to have an actuator with a fast response so that the dynamics of the actively 
controlled airplane will be the dominant component of the system response. A slow actuator is akin to a time-delay that can cause 
performance and stability problems. 

The second-order system approximation has allowed us to gain insight into the relationship between the parameter K and the 
system response, as measured by percent overshoot and time-to-peak. Of course, the gain K — 0.16 is only a starting point in the 
design. Why? Because we in fact have a third-order system and must consider the effect of the third pole (which we have ignored 
so far). 

We can develop a simulation of the control system to quantify the performance. We can readily vary important system parameters 
(such as K) and check the resulting performance (as measured by step response characteristics). LabVIEW provides a graphical 
methodology for developing the simulation using block diagrams. As shown in Fig. 4.8, LabVIEW presents the window of the 
system in block diagram form. 


l Aircraft Control and Simulation (2nd Edition) by Brian L. Stevens and Frank L. Lewis, Wiley-Interscience, 2003. 
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Figure 4.8: LabVIEW simulation window with aircraft, controller, and actuator dynamics. 


Before starting the simulation we can opt to change many of the parameters of the system we are modeling. To change any part 
in the model, simply open (double-click with the mouse button) the Controller Gain, Aileron Actuator, Aircraft Dynamics, or Gyro 
Gain and enter new parameters. 


For example, opening the Controller Gain displays the window shown in Fig. 4.9. We use this window to specify the value of the 
gain introduced into the system. To move the window and avoid obstructing the diagram, we can click and hold the mouse button 
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Figure 4.9: Pop-up window used to change the controller gain. 
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Figure 4.10: Pop-up window used to change the aileron actuator transfer function. 


anywhere on the title bar and we drag the window some place else. As we change the parameters they are immediately updated on 
the diagram. Figures 4.10 and 4.11 show the display windows opened for the aileron actuator and the aircraft dynamics. 

We need to set a few more parameters before the model is ready for simulation. Figure 4.12 shows the Configure Simulation 
Parameters window that allows users to select different numerical integration routines and the parameters for the initial and final 
time, maximum and minimum step size, and relative and absolute tolerances. 

After making the desired changes, select the OK button. During the simulation execution, the time response plot is dynamically 
updated. The time response should now resemble Fig. 4.13. 

To fine-tune the simulation, we can change any or all of the transfer functions or gains by the steps that we just covered. At 
this point the LabVIEW setup has K — 0.16. From our previous analytic analysis, we determined that as K decreased, the percent 
overshoot also decreased, while the time-to-peak simultaneously increased. With an interactive simulation we can verify this trend, 
except that now we have the original third-order system in the simulation rather than the approximate second-order system on which 
the analytic analysis was based. With the second-order system approximation, we estimate that with K = 0.16 the percent overshoot, 
P.O., is 20% and the time-to-peak, Tp, is 2.62 seconds. What is the system performance with the third-order system? LabVIEW 
provides a way to investigate this easily. The results of the simulation are shown in Fig. 4.13. The percent overshoot is slightly more 
than 2046 due to the presence of a third pole that we ignored in the second-order system analytic analysis. In fact, with K — 0.16 
we have 

P.O. = 20.590 and T, = 2.73 seconds. 


The percent overshoot is reduced from 20.5% (with K = 0.16) to 9.5% (with K = 0.1). Also the time-to-peak increases from 2.73 
seconds (with K — 0.16) to 3.74 seconds (with K — 0.1). We can verify that when K — 0.2, the percent overshoot and time-to-peak 
are 26.5% and 2.38 seconds, respectively. © 
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Figure 4.11: Pop-up window used to change the aircraft dynamics transfer function. 
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Figure 4.12: Selecting Parameters from the Configure Simulation Paraeters menu. 
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Figure 4.13: The LabVIEW simulation with the scaled graph showing the aircraft bank angle response to a step input. 


Stability of Linear Feedback Systems 


This chapter begins with a discussion of the Routh-Hurwitz stability method. We will see how LabVIEW can assist us in stability 
analysis by providing a readily accessible method for computing the poles of the characteristic equation. For the case when the 
characteristic equation is a function of a single parameter, it will be possible to generate a plot displaying the movement of the poles 
as the parameter varies. The chapter concludes with discussion on stability for state variable systems. 


5.1 Routh—Hurwitz Stability 


The Routh—Hurwitz criterion is a necessary and sufficient criterion for stability. Given a characteristic equation with fixed coefficients, 
we can use Routh—Hurwitz to determine the number of roots in the right half-plane. For example, consider the characteristic equation 


g(s) = S +57 2s 4-24 20 


associated with the closed-loop control system shown in Fig. 5.1. The corresponding Routh array is shown in Fig. 5.2. The two 
sign changes in the first column indicate that there are two roots of the characteristic polynomial in the right half-plane; hence 
the closed-loop system is unstable. Using LabVIEW, we can verify the Routh-Hurwitz result by directly computing the roots of 
the characteristic equation using the CD Poles function (as shown in Fig. 5.3). Recall that the CD Poles function computes the 
system poles. Whenever the characteristic equation is a function of a single parameter, the Routh-Hurwitz method can be utilized 
to determine the range of values that the parameter may take while maintaining stability. Consider the closed-loop feedback system 
in Fig. 5.4. The characteristic equation is 


g(s) = s? +2s* c As + K — 0. 


I 


G(s) — —3 Y(s) 


s? * s? * 25 + 23 





Ist sign change 






2nd sign change 


Figure 5.2: Routh array for the closed-loop control system with T (s) = Y (s)/ R(s) = 1/ (s? +57 + 2s + 24). 


39 


40 CHAPTER 5 STABILITY OF LINEAR FEEDBACK SYSTEMS 





Figure 5.3: Using the CD Poles function to compute the closed-loop control system poles of the system shown in Fig. 5.1. 


The associated Routh array is given by 


S 
S 


N 
N e 
O rN A 


S y 
1| K 


The Routh-Hurwitz criterion states that for closed-loop stability, all terms in the first column of the Routh array must be positive 
Thus we have 


8$- K>0 => K <8, 
K>0 => K -0. 


Y(s) 





Figure 5.4: Closed-loop control system with T (s) = Y (s)/ R(s) = 1/ (s? + 25? -- As + K). 
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Figure 5.5: Plot of root locations of g(s) = s? +28? +45 + K for0 < K < 20. 


We can use LabVIEW to verify this result graphically. As shown in Fig. 5.5, we establish a vector of values for K at which we 
wish to compute the roots of the characteristic equation. Then using the Polynomial Roots function, we calculate and plot the roots 
of the characteristic equation, as shown in Fig. 5.5. It can be seen that as K increases, the roots of the characteristic equation move 
toward the right half-plane as the gain tends toward K = 8, and eventually into the right half-plane when K > 8. 


The VI in Fig. 5.6 uses a For Loop function. The For Loop structure is described in more detail in Fig. 5.6. This structure 
provides a mechanism for repeatedly executing a series of statements a given number of times. The For Loop structure sets up a 
repeating calculation loop. Figure 5.6 describes the For Loop structure format and provides an illustrative example of its usefulness. 
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Figure 5.6: Using the For Loop structure. 


The example sets up a loop that repeats ten times. During the ith iteration, where 1 < i < 10, the ith element of the vector a is set 
equal to 20, and the scalar b is recomputed. 

The Routh-Hurwitz method allows us to make definitive statements regarding absolute stability of a linear system. The method 
does not address the issue of relative stability, which is directly related to the location of the roots of the characteristic equation. 
Routh-Hurwitz tells us how many poles lie in the right half-plane, but not the specific location of the poles. With LabVIEW, we can 
easily calculate the poles explicitly, thus allowing us to investigate relative stability. 


Example 5.1 Tracked Vehicle Control 


The design objective is to find a and K such that the system is stable and the steady-state error for a ramp input is less than or equal to 
24% of the command. We can use the Routh-Hurwitz method to aid in the search for appropriate values of a and K. The closed-loop 
characteristic equation is 


q(s) = sf 4-8s? +1757 + (K +10)s +aK — 0. 
Using the Routh array, we find that for stability we require 


K « 126, aK > 0. 
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Figure 5.7: Stability region for a and K for two-track vehicle turning control. 


For positive K it follows that we can restrict our search to 0 < K < 126 and a > 0. Our approach will be to use LabVIEW to 
find a parameterized a versus K region in which stability is assured. Then we can find a set of (a, K) belonging to the stable region 
such that the steady-state error specification is met. This procedure, shown in Fig. 5.7, involves selecting a range of values for a 
and K and computing the roots of the characteristic polynomial for specific values of a and K. For each value of K, we find the 
first value of a that results in at least one root of the characteristic equation in the right half-plane. The process is repeated until the 
entire selected range of a and K is exhausted. The plot of the (a, K) pairs defines the separation between the stable and unstable 
regions. The region to the left of the plot of a versus K in Fig. 5.7 is the stable region. If we assume that r(t) = At, t > O, then the 
steady-state error is 


s(s + 1)(s + 2)(s +5) A 10A 
ess = lim s 2 __. i 
s>0 s(st1)(s+2)(s+5)+K(sta) s aK 


where we have used the fact that 


OE NN as, 
1+ G.G(s) s(s + Ds +2)(s +5) + K(s +a) 
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Figure 5.8: Ramp response for a = 0.6 and K = 70 for two-track vehicle turning control. 


Given the steady-state specification, ess < 0.24A, we find that the specification is satisfied when 


10A 
— < 0.24A, 
aK 


or 


aK > 41.67. (5.1) 


Any values of a and K that lie in the stable region in Fig. 5.7 and satisfy Eq. (5.1) will lead to an acceptable design. For example, 
K = 70 and a = 0.6 will satisfy all the design requirements. The closed-loop transfer function (with a = 0.6 and K = 70) is 


70s + 42 


T(s) = —— SS —.. 
(s) s^ + 8s? + 17s? + 80s + 42 


The associated closed-loop poles are 


=  —7.0767, 

=  —0.5781, 

—0.1726 + 3.19957, and 
=  —0.1726 — 3.19957. 


Annn 
| 


The corresponding unit ramp input response is shown in Fig. 5.8. The steady-state error is less than 0.24, as desired. o 
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Figure 5.9: Computing the characteristic polynomial of A with the CD Characteristic Polynomial function. 


5.2 The Stability of State Variable Systems 


Now let us turn to determining the stability of systems described in state variable form. Suppose we have a system in state-space 
form 
x = Ax+ Bu. 


The stability of the system can be evaluated with the characteristic equation associated with the system matrix A. The characteristic 
equation is 
det(sI — A) = 0. (5.2) 


The left-hand side of the characteristic equation is a polynomial in s. If all of the roots of the characteristic equation have negative 
real parts (1.e., Re(s;) < 0), then the system is stable. 

When the system model is given in state variable form, we must calculate the characteristic polynomial associated with the A 
matrix. In this regard we have several options. We can calculate the characteristic equation directly from Eq. (5.2) by manually 
computing the determinant of sI — A. Then we can compute the roots using the polynomial roots function to check for stability, or 
alternatively, we can utilize the Routh-Hurwitz method to detect any unstable roots. Unfortunately the manual computations can 
become lengthy, especially if the dimension of A is large. We would like to avoid this manual computation if possible. As it turns 
out, LabVIEW can assist in this endeavor. The Matrix Characteristic Polynomial function is used to compute the characteristic 
equation associated with A, as illustrated in Fig. 5.9. 

The input matrix A is 


—8 -16 -6 
A= 1 0 0|, 
0 I 0 


and the associated characteristic polynomial is 
s? 4- 8s? + 16s -- 6 — 0. 
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If A is an n x n matrix, Matrix Characteristic Polynomial(A) is the characteristic equation represented by the n + 1 element row 
vector whose elements are the coefficients of the characteristic equation. 


Example 5.2 Stability Region for an Unstable Plant 


A jump-jet aircraft has a control system as shown in Fig. 5.10. Assume that z > 0 and p > 0. The system is open-loop unstable 
(without feedback) since the characteristic equation of the plant and controller is 


s(s — D) (s + p) = s[s? + (p — Ds — p] = 0. 


Note that since one term within the bracket has a negative coefficient, the characteristic equation has at least one root in the right-hand 
s-plane. The characteristic equation of the closed-loop system is 


s? -(p— Ds? + (K — p)s + Kz — 0. 


The goal is to determine the region of stability for K, p, and z. The Routh array is 


S l K — p 
s | p—1 Kz 
S bz 0 

1 Kz 


where 
_ (p - Y(K = p) - Kz 


b 
2 PN 


From the Routh-Hurwitz criterion, we require Kz > 0 and p > 1. Setting b» = 0, we have 


(p—-1)(K — p)-Kz=K [(p-1)-z]-—p(p—-) =0. 


Therefore we require that 


Ex p(p — 1) (5.3) 


(p-1-z 
Consider three cases: 
l. z > p—1: Since p > 1, any K > O will satisfy the stability conditions. 


2. z= p — 1: Thereisno0 < K < oo that leads to stability. 


3. z< p—1: AnyO « K < œ satisfying stability condition Eq. (5.3) for a given p and z will result in stability. 





















K(s + 1 
R(s) see) — Y) 
Aircraft Bom eet) Actual 
desired orientation 


orientation 


Figure 5.10: Control system for jump-jet aircraft. Assume that z > 0 and p > Q. 
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Figure 5.11: The three-dimensional region of stability lies above the surface shown. 


The stability conditions can be depicted graphically. The LabVIEW script used to generate a three-dimensional stability surface 
is shown in Fig. 5.11. One acceptable stability point is z — 1, p — 10, and K — 15. © 


Root Locus Method 


LabVIEW can readily be used to obtain accurate root locus plots. The fundamental concepts associated with sketching a root locus 
are embedded in the CD Root Locus function. If you are not familiar with how to sketch a root locus, it is recommended that you 
learn. It is essential for good control system design with root locus to fully understand the underlying theory. The chapter begins 
with a discussion on obtaining a root locus plot with LabVIEW. This is followed by a discussion of the connections between the 
partial fraction expansion, dominant poles, and the closed-loop system response. Root sensitivity is covered in the final paragraphs. 
The functions covered in this section are CD Root Locus and Partial Fraction Expansion. 


6.1 Obtaining a Root Locus Plot 
Consider the closed-loop control system in Fig. 6.1. The closed-loop transfer function is 


OYW K+ 6 +3) 
© R(s)) ss +2)(8 +3) 4+ K(s +1) 





T (s) 


The characteristic equation can be written as 
s+l 
K —————— = 
s(s + 2)(s + 3) 
The characteristic equation in Eq. (6.1) is in the necessary Evans form for generating root locus plots using LabVIEW. The general 
form of the characteristic equation necessary for application of the CD Root Locus function is 


(6.1) 


p(s) — 
q(s) 
where K is the parameter of interest to be varied from 0 < K « oo. The CD Root Locus function is used to obtain the root locus 


plot associated with Eq. (6.1). For the characteristic equation in Eq. (6.1), the associated root locus plot is shown in Fig. 6.2. 
The steps required to use CD Root Locus to obtain a root locus plot with LabVIEW are as follows: 


| c-K 0, (6.2) 


1. Obtain the characteristic equation in the form given in Eq. (6.2), where K is the parameter of interest. 
2. Use the CD Root Locus function to generate the plots. 


Referring to Fig. 6.2, we see that as K increases, two branches of the root locus break away from the real axis. This means that 
for some values of K the closed-loop system characteristic equation will have complex roots. For example, we find that when 


K(s +1) 
s(s +2) 


1 
S 3 











Y(s) 





Figure 6.1: Closed-loop system. 
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Figure 6.2: The CD Root Locus function. 


K = 20.5775, the closed-loop transfer function has three poles and two zeros at poles: 


—2.0505 + 4.32271 


poles: s = | —2.0505 — 4.32271 , zeros: s = ( | ) ; 
—0.8989 
Considering the closed-loop pole locations only, we would expect that the real pole at s = —0.8989 would be the dominant pole. 
To verify this, we can study the closed-loop system response to a step input, R(s) = 1/s. For a step input we have 
20.5775 1 3 1 
Y(s) = WSs t D(st3) — S (6.3) 


s(s -- 2)(s +3) +20.5775(s +1) s 


Generally the first step in computing y(t) is to expand Eq. (6.3) in a partial fraction expansion. The Partial Fraction Expansion 
function can be used to expand Eq. (6.3) in a partial fraction expansion, as shown in Fig. 6.3. The partial fraction expansion of 
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Figure 6.3: Partial fraction expansion of Eq. (6.3). 


Eq. (6.3) 1s 
— ].3786 + 1.70107 — 1.3786 — 1.70107 —0.2429 3 


¥(s) = — A 5 
s + 2.0505 + 4.32281 s + 2.0505 — 4.3228: s + 0.8989 


Comparing the terms of the partial fraction expansion, we see that the coefficient of the term corresponding to the pole at s = —0.8989 
is considerably smaller than the coefficient of the terms corresponding to the complex-conjugate poles at s = —2.0505 + j4.3227. 
From this we expect that the influence of the pole at s = —0.8989 on the output response y(t) does not dominate. The settling time 
to within 2% of the final value is then predicted by considering the complex-conjugate poles. The poles at s = —2.0505 + j4.3227 
correspond to a damping of ¢ = 0.4286 and a natural frequency of œn = 4.7844. Thus the settling time is predicted to be 


= 1.95 seconds. 


N 
Q 





n 


Using the CD Paraetric Time Response function, as shown in Fig. 6.4, we find that T, ~ 1.6 seconds. Hence our approximation 
of settling time T, ~ 1.95 is a fairly good approximation. 

In this example the role of the system zeros on the transient response is illustrated. The proximity of the zero at s — —1 to the 
pole at s = —0.8989 reduces the impact of that pole on the transient response. The main contributors to the transient response are 
the complex-conjugate poles at s = —2.0505 + 74.3228 and the zero at s = —3. 
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Figure 6.4: CD Step Response function with the closed-loop system in Fig. 6.1 with K = 20.5775. 


6.1.1 Sensitivity and the Root Locus 


The roots of the characteristic equation play an important role in defining the closed-loop system transient response. The effect of 
parameter variations on the roots of the characteristic equation is a useful measure of sensitivity. The root sensitivity can be defined 


to be 
li Or; 


K JK/K` (om) 


S 
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Figure 6.5: Sensitivity calculations for the root locus for a 596 change in K. 


We can utilize Eq. (6.4) to investigate the sensitivity of the roots of the characteristic equation to variations in the parameter K. If 
we change K by a small finite amount AK, and evaluate the modified root r; + Ari, it follows that 


sa ON (6.5) 
K” AK/K’ ' 


In general, the quantity Se is a complex number. Referring back to the third-order example of Fig. 6.1, if we change K by a factor 
of 5%, we find that the dominant complex-conjugate pole at s = —2.0505 + j4.3228 changes by 


Ar; = —0.0025 — 0.11681 


when K changes from K = 20.5775 to K = 21.6064. From Eq. (6.5), it follows that 


, . —0.0025 — 0.1168: | 
S m — — — —0.0494 — 2.3355i. 
1.0289/20.5775 


The sensitivity S can also be written in the form 
S = 2.3360 /268.7872". 


The magnitude and direction of S e provides a measure of the root sensitivity. The script used to perform these sensitivity calculations 
is shown in Fig. 6.5. The root sensitivity measure is useful for comparing the sensitivity for various system parameters at different 
root locations. 


CHAPTER 7 


Frequency Response Methods 


This chapter begins with an introduction to the Bode plot and then discusses the connection between the frequency response and 
performance specifications in the time domain. We conclude with an illustrative example of designing a control system in the 
frequency domain. The LabVIEW function covered is CD Bode. The CD Bode function is used to generate a Bode plot. We will 
use transfer function models in this chapter, but the use of the CD Bode function is exactly the same for state variable models as 
with transfer functions, except that the input is a state-space object instead of a transfer function object. 


7.1 Bode Plots 


Consider the transfer function 


S(1 + 0.15) 


G(s) = ——__________.. 
s(1+0.5s)(1 + Ss + <552) 


(7.1) 


The Bode plot corresponding to Eq. (7.1) is shown in Fig. 7.1. The plot consists of the logarithmic gain in dB versus o in one plot 
and the phase $ (c) versus w in a second plot. As with the root locus plots, it will be tempting to rely exclusively on LabVIEW to 
obtain your Bode plots. Treat LabVIEW as one tool in a tool kit that can be used to design and analyze control systems. It is essential 
to develop the capability to sketch Bode plots. There is no substitute for a clear understanding of the underlying theory. 


Keeping in mind the goal of designing control systems that satisfy certain performance specifications given in the time domain, 
we must establish a connection between the frequency response and the transient time response of a system. The relationship between 
specifications given in the time domain to those given in the frequency domain depends upon approximation of the system by a 
second-order system with the poles being the system dominant poles. 


Consider the second-order system shown in Fig. 7.2. The closed-loop transfer function is 


w2 


T(s) = ———— ———. Ju 
‘s) s? + 2t£cons + o? ae) 


The Bode plot magnitude associated with the closed-loop transfer function in Eq. (7.2) is shown in Fig. 7.3. The relationship between 
the resonant frequency, w, the maximum of the frequency response, M p„, and the damping ratio, ¢, and the natural frequency, cy, 
is shown in Fig. 7.4. The information in Fig. 7.4 will be quite helpful in designing control systems in the frequency domain while 
satisfying time-domain specifications. The relationship between the resonant frequency, œw, and the natural frequency, œn, is 


Wy -—41—2£?o0,, € «0.707 (7.3) 
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Figure 7.1: The Bode plot associated with Eq. (7.1). 
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Figure 7.2: A second-order closed-loop system. 
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Figure 7.3: Magnitude characteristic of the second-order system. 
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Figure 7.5: (a) Engraving machine control system. (b) Block diagram model. 


for second-order systems. The relationship between the maximum frequency response, M y, , and the damping ratio, ¢, is 


My, = (2 1 = g) | £ «0707 (7.4) 


Example 7.1 Engraving Machine System 


Consider the block diagram model in Fig. 7.5. Our objective is to select K so that the closed-loop system has an acceptable time 
response to a step command. A functional block diagram describing the frequency-domain design process is shown in Fig. 7.6. First 
we choose K — 2 and then iterate on K if the performance is unacceptable. The script shown in Fig. 7.7 is used in the design. The 
value of K is defined at the command level. Then the script is executed and the closed-loop Bode diagram is generated. The values 
of Mp, and œ, are determined by inspection from the Bode plot. Those values are used in conjunction with Fig. 7.4 to determine the 
corresponding values of ¢ and œn. Given the damping ratio, ¢, and the natural frequency, wn, the settling time and percent overshoot 
are estimated using the formulas 


4 — N 
T, © , P.O. 100exp 


COn J/1- zz 
If the time-domain specifications are not satisfied, then we adjust K and iterate. The values for ¢ and cj, corresponding to K = 2 are 


¢ = 0.29 and œn = 0.88. This leads to a prediction of P.O. = 37% and T; = 15.7 seconds. The step response, shown in Fig. 7.8, 
is a verification that the performance predictions are quite accurate and that the closed-loop system performs adequately. © 





In this example, the second-order system approximation is reasonable and leads to an acceptable design. However, the second- 
order approximation may not always lead directly to a good design. Fortunately with LabVIEW we can construct an interactive 
design facility to assist in the design process by reducing the manual computational loads while providing easy access to a host of 
classical and modern control tools. 
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Figure 7.6: Frequency design functional block diagram for the engraving machine. 
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Figure 7.7: Script for the design of an engraving machine. 
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Figure 7.8: Engraving machine step response for K — 2. 


Stability in the Frequency Domain 


This chapter considers the Nyquist diagram and the Nichols chart. Two examples are used illustrate the frequency-domain design 
approach. We also present an illustrative example that shows how to deal with a time delay in the system utilizing a Padé approximation. 
The LabVIEW functions covered in this chapter are CD Nyquist, CD Nichols, CD Gain and Phase Margin, and CD Construct 
Special TF Model. 


8.1 Nyquist Plots 


It is generally more difficult to generate the Nyquist plot manually than the Bode plot. However, we can use LabVIEW to obtain the 
Nyquist plot. The Nyquist plot is generated with the CD Nyquist function, as shown in Fig. 8.1. 
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Figure 8.1: The CD Nyquist function. 
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Figure 8.2: A closed-loop control system example for relative stability analysis. 


Relative stability measures of gain margin and phase margin can be determined from both the Nyquist plot and the Bode plot. 
The gain margin is a measure of how much the system gain would have to be increased for the loop transfer function L(j@) locus 
to pass through the (—1, 0) point, thus resulting in an unstable system. The phase margin is a measure of the additional phase lag 
required before the system becomes unstable. Gain and phase margins can be determined from both the Nyquist plot and the Bode 


plot. 


Consider the system shown in Fig. 8.2. Relative stability can be determined from the Bode diagram using the CD Gain and 
Phase Margin function, which is shown in Fig. 8.3. The VI used to generate the Nyquist plot for the system in Fig. 8.2 1s shown 
in Fig. 8.4. In this case, the number of poles of the loop transfer function L(s) with positive real parts is zero, and the number of 
counterclockwise encirclements of —1 is zero; hence the closed-loop system is stable. We can also determine the gain margin and 


phase margin, as indicated in Fig. 8.4. 
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Figure 8.3: The Bode diagram for the system in Fig. 8.2 with the gain margin and phase margin. 
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Figure 8.4: The Nyquist plot for the system in Fig. 8.2 with gain and phase margins. 


$.2 Nichols Charts 


Nichols charts can be generated using the CD Nichols function, shown in Fig. 8.5. The Nichols chart shown in Fig. 8.6 is for the 
system 
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Figure 8.5: The CD Nichols function. 
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Figure 8.6: Nichols chart for the system in Eq. (8.1). 
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Figure 8.7: Liquid level control system block diagram. 


Example 8.1 Liquid Level Control System 


Consider a liquid level control system described by the block diagram shown in Fig. 8.7. Notice that this system has a time delay. 
The loop transfer function is given by 


31.5 Ls 


L(s) = ———_______.——_ exp’ . (8.2) 
(s + 1)30s+ DOE +5341) 

We first rearrange Eq. (8.2) in such a way that L (5) has a transfer function form with polynomials in the numerator and denominator. 

To do this we can make an approximation to e~*/ with the CD Construct Special TF Model function. The CD Construct Special 

TF Model function is shown in Fig. 8.8 with the time delay as an input. For example, suppose our time delay is T = 1 second and 

we want a second-order approximation n — 2. Then, using the CD Construct Special TF Model function we find that 


_ 0.074357 — 0.44605 + 0.8920 
0.07435? + 0.4460s + 0.8920" 


—s 


(8.3) 


Substituting Eq. (8.3) into Eq. (8.2), we have 


31.5(0.0743s? — 0.4460s + 0.8920) 


L(s) = —-?ASN , 
(s + 1)(30s + D qw 1)(0.0743s? + 0.44605 + 0.8920) 

Now we can build a script to investigate the relative stability of the system using the Bode diagram. Our goal is to have a phase 

margin of 30°. The associated script is shown in Fig. 8.8. To make the VI interactive, we let the gain K (now set at K = 31.5) 

be adjustable. We set K and run the VI to check the phase margin and iterate, if necessary. The final selected gain is K = 16. 

Remember that we have utilized a second-order Padé approximation of the time delay in our analysis. Q 
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num = K 

denl = [1 1]; 

den? = [301]; 

den3 = [1/9 1/3 1]; 
sysl = tf(num, denl); 
sys? = tf(1,den2); 
sys3 = tf(1,den3); 

G = sysl*sys2*sys3 
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Figure 8.8: Bode diagram for the liquid level control system. 
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Figure 8.9: Remotely controlled reconnaissance vehicle's speed control system. 


Example 8.2 Remotely Controlled Reconnaissance Vehicle 


Consider the speed control system for a remotely controlled reconnaissance vehicle shown in Fig. 8.9. The design objective is to 
achieve good control with low steady-state error and low overshoot to a step command. Building a VI will allow us to perform many 
design iterations quickly and efficiently. First we investigate the steady-state error specification. The steady-state error, e;;, to a unit 
step command is 


1 


ess 
The effect of the gain K on the steady-state error is clear from Eq. (8.4): If K = 20, the error is 9% of the input magnitude; if 


K = 10, the error is 17% of the input magnitude. Now we can investigate the overshoot specification in the frequency domain. 
Suppose we require that the percent overshoot be less than 50%. Solving 


P.O. © 100 exp 57/ V 1-€ < 50 


for £ yields ¢ > 0.215. For second-order systems we have the relationship between M ,, and ¢ given by 


My, = (2c VA E gy | (8.5) 


With ¢ > 0.215, we find that Mp, < 2.45. We must keep in mind that the formula in Eq. (8.5) is exact for second-order systems 
only and can be used here only as a guideline. We now compute the closed-loop Bode diagram and check the values of Mp,- Any 
gain K for which M,, < 2.45 may be a valid gain for our design, but we will have to investigate further to include step responses 
to check the actual overshoot. The script in Fig. 8.10 aids us in this task. We can investigate further the gains K — 20, 10, and 4.44 
(even though M,, > 2.45 for K = 20). We plot the step responses to quantify the overshoot as shown in Fig. 8.11. Additionally, 
we could have used a Nichols chart to aid the design process as shown in Fig. 8.12. The results of the analysis are summarized in 
Table 8.1 for K = 20, 10, and 4.44. 


We choose K = 10 as our design gain. The gain margin is GM = 49.56 dB and the phase margin is PM = 26.11°. © 
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Table 8.1: Compensator Design Results 


K 4.44 10 20 

Percent overshoot 32.44 48.4 61.4 
Settling time (seconds) 4.94 5.46 6.58 
Peak time (seconds) 1.19 0.88 0.67 
E 31% 16.7% 9.1% 
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Figure 8.10: Remotely controlled vehicle closed-loop system Bode diagram 
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Figure 8.11: Remotely controlled vehicle step response. 
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Figure 8.12: Remotely controlled vehicle Nichols chart. 


Design of Feedback Control Systems 


In this chapter the compensation of control systems is illustrated using frequency response and root locus methods. We reconsider 
a rotor winder design example to illustrate the use of LabVIEW scripts in designing and developing control systems with good 
performance characteristics. Lead and lag compensators are examined for this design example and the system response is obtained 
using LabVIEW. 


Example 9.1 Rotor Winder Control System 


Consider the rotor winder control system shown in Fig. 9.1. The design objective is to achieve high steady-state accuracy to a ramp 
input. The steady-state error to a unit ramp input, R(s) = 1/s?, is 


1 
e =; 
O K, 


where 

" Gc(s) . Ge(s) 

im $————— —— —— = lim 
s0 s(s+5)(s+10) sso 50 
The performance specification of overshoot and settling time must be considered as well as steady-state tracking error. In all 
likelihood, a simple gain will not be satisfactory, so we will also consider controllers utilizing lead and lag compensators using both 
Bode diagram and root locus plot design methods. Our approach is to develop a series of LabVIEW scripts to aid in the compensator 
designs. Consider first a simple gain controller, G-(s), where 


Ky = 





G.e(s) = K. 
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Figure 9.1: (a) Rotor winder control system. (b) Block diagram. 
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Figure 9.2: Transient response for simple gain controller. 


Then, the steady-state error is 
50 


ess = K` 
The larger we make K, the smaller is the steady-state error ess. However, we must also consider the effect that increasing K has on 
the transient response, as shown in Fig. 9.2. When K = 500, our steady-state error for a ramp is 10%, but the overshoot is 70% and 
the settling time is approximately 8 seconds for a step input. This is unacceptable performance and thus we turn to compensation. 
The two important compensator types considered are lead and lag compensators. First, consider the lead compensator 


K(s +z) 


Ge(s) = (s+ p) , 


where |z| « |p|. The lead compensator provides the capability to improve the transient response. We will use a frequency-domain 
approach to design the lead compensator. We desire a steady-state error of less than 10% to a ramp input, or K, = 10. In addition to 
the steady-state specifications, we desire to meet certain performance specifications: (1) settling time (2% criterion) T; < 3 seconds, 
and (2) percent overshoot for a step input P.O. «10946. Solving for ¢ and œn using 


P.O. =100expS7/V!-* —10 and T= Eu 





T2 
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yields ¢ = 0.59 and @, = 2.26. We thus obtain the phase margin requirement: 


C 
x — — 60°. 
P pm 0.01 


The steps leading to the final design are as follows: 


L. 


Zz; 


ie 


8. 


Obtain the uncompensated system Bode diagram with K = 500, and compute the phase margin. 
Determine the amount of necessary phase lead. 
Evaluate o from sin $,, = (a — 1)/(o + 1). 


Compute 10 log œ and find the frequency wn on the uncompensated Bode diagram where the magnitude curve is equal to -10 
log o. 


In the neighborhood of w, on the uncompensated Bode, draw a line through the 0-dB point at œn with slope equal to the current 
slope plus 20 dB/decade. Locate the intersection of the line with the uncompensated Bode to determine the lead compensation 
zero location. Then calculate the lead compensator pole location as p — oz. 


. Obtain the compensated Bode and check the phase margin. Repeat any steps if necessary. 


Raise the gain to account for attenuation (1/o). 


Verify the final design with simulation using step functions, and repeat any steps if necessary. 


Three scripts are utilized in the design. The design scripts are shown in Figs. 9.3-9.5. The first script is for the Bode plot of the 
uncompensated system, as shown in Fig. 9.3. The detailed Bode plot of the compensated system is shown in Fig. 9.4. Figure 9.5 
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Figure 9.3: Bode diagram. 
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Figure 9.4: Lead compensator Bode diagram. 
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Figure 9.5: Lead compensator step response. 
displays the LabVIEW script for the step response analysis. The final lead compensator design is 
G(s) 1800(s + 3.5) 
s) = ———————, 
É (s + 25) 
where K = 1800 was selected after iteratively using the LabVIEW script. The settling time and overshoot specifications are 


satisfied, but K, = 5, resulting in a 20% steady-state error to a ramp input. It is possible to continue the design iteration and refine 
the compensator somewhat, although it should be clear that the lead compensator has added phase margin and improved the transient 


response as anticipated. 


To reduce the steady-state error, we can consider the lag compensator, which has the form 


K(s +z) 


G.(s) = (s+ p) , 


where |p| « |z|. We will use a root locus approach to design the lag compensator, although it can be done using Bode as well. The 
desired root location region of the dominant roots is specified by ¢ = 0.59 and œn = 2.26. The steps in the design are as follows: 


1. Obtain the root locus of the uncompensated system. 


2. Locate suitable root locations on the uncompensated system that lie in the region defined by ¢ = 0.59 and w, = 2.26. 


3. Calculate the loop gain at the desired root location and the system error constant, K 


4. Compute a = K 


Uuncomp * 


Vcomp / K vuncomp where K veomp = 10. 
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5. With a known, determine suitable locations of the compensator pole and zero so that the compensated root locus still passes 
through desired location. 


6. Verify with simulation and repeat any steps if necessary. 


The design methodology is illustrated in Figs. 9.6-9.8. We can compute the gain K associated with the roots of our choice on the 
uncompensated root locus that lie in the performance region. We then compute o to ensure that we achieve the desired K,. We 
place the lag compensator pole and zero to avoid impacting the uncompensated root locus. In Fig. 9.7, the lag compensator pole and 
zero are very near the origin, at z = —0.1 and p = —0.01. The settling time and overshoot specifications are nearly satisfied, and 
K, = 10, as desired. It is possible to continue the design iteration and refine the compensator somewhat, although it should be clear 
that the lag compensator has improved the steady-state errors to a ramp input relative to the lead compensator design. The final lag 
compensator design is 
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Figure 9.6: Lag compensator: uncompensated root locus. 
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Figure 9.7: Lag compensator: compensated root locus. 
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Figure 9.8: Lag compensator: Step response. 
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Design of State Variable Feedback Systems 


Controllability and observability of a system in state variable feedback form can be checked using the LabVIEW functions CD 
Controllability Matrix and CD Observability Matrix, respectively. The inputs to the CD Controllability Matrix function are the 
system matrix A and the input matrix B; the output of CD Controllability Matrix is the controllability matrix Pe. Similarly, the 
input to the CD Observability Matrix function is the system matrix A and the output matrix C and the output of CD Observability 
Matrix is the observability matrix P,. Note that the controllability matrix, P,, is a function only of A and B, while the observability 
matrix, Po, is a function only of A and C. 


Example 10.1 Satellite Trajectory Control 


Consider a satellite in a circular, equatorial orbit at an altitude of 250 nautical miles above the earth, as illustrated in Fig. 10.1. The 
satellite motion (in the orbit plane) is described by the normalized state variable model 


0 1 0 0 0 0 

, 307 0 0 2o 1 0 

X es 0 0 0 1 x + 0 uy + o |": (10.1) 
0 -20 0 0 0 1 


where the state vector x represents normalized perturbations from the circular, equatorial orbit, u, is the input from a radial thruster, 
ur is the input from a tangential thruster, and œ = 0.0011 rad/s (approximately one orbit of 90 minutes) is the orbital rate for the 


N 





Circular, equitorial orbit 


Figure 10.1: The satellite in an equatorial, circular orbit. 
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Figure 10.2: Controllability with radial thrusters only. 


satellite at the specific altitude. In the absence of perturbations, the satellite will remain in the nominal circular, equatorial orbit. 
However, disturbances such as aerodynamic drag can cause the satellite to deviate from its nominal path. The problem is to design a 
controller that commands the satellite thrusters in such a manner that the actual orbit remains near the desired circular orbit. Before 
commencing with the design, we check controllability. In this case we investigate controllability using the radial and tangential 
thrusters independently. Suppose the tangential thruster fails (1.e., ur = 0), and only the radial thruster is operational. Is the satellite 
controllable from u, only? We answer this question by using LabVIEW to determine the controllability. Using the script shown in 
Fig. 10.2, we find that the determinant P, is zero; thus the satellite is not completely controllable when the tangential thruster fails. 


Suppose now that the radial thruster fails (1.e., u, = 0) and that the tangential thruster is functioning properly. Is the satellite 
controllable from u; only? Using the script in Fig. 10.3, we find that the satellite is completely controllable using the tangential 
thruster only. o 


We conclude this section with a controller design for an automatic test system using state variable models. The design approach 
utilizes root locus methods and incorporates LabVIEW scripts to assist in the procedure. 
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Figure 10.3: Controllability with tangential thrusters only. 


Example 10.2 Automatic Test System 


The state-space representation for an automatic test system is 


X = AX + Bu, (10.2) 
where 
0 1 0 0 
A—| 0 -1 1 |, B = 0 
0 0 —5 K 


Our design specifications are a step response with (1) a settling time (2% criterion) less than 7, < 2 seconds and (2) an overshoot 
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less than P.O. < 4%. We assume that the state variables are available for feedback, so that the control is given by 
u = (—K| -K — K3) x = Kx. (10.3) 
We must select the gains K, K1, K2, and K3 to meet the performance specifications. Using the design approximations 


4 


Cn 





P.O. = 100exp *"/v 1—6? <4 and = ae 


we find that 
¢ > 0.72 and Wy, > 2.8. 


This defines a region in the complex plane in which our dominant roots must lie, so that we expect to meet the design specifications, 
as shown in Fig. 10.4. Substituting Eq. (10.3) into Eq. (10.2) yields 


0 ] 0 
xm 0 —] 1 x = Hx, (10.4) 
—KK, -KK —(5+KK3) 
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Figure 10.4: Root locus for the automatic test system. 
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where H = A — BK. The characteristic equation associated with Eq. (10.4) can be obtained by evaluating det(sI — H) =0, resulting 
in 


(10.5) 


K3+K K 
s(s + D(s +5) 4+ K Ka (^. SE, J — 0. 


K3 K3 
If we view K K3 as a parameter and let Kı = 1, then we can write Eq. (10.5) as 
s? -F MESS T 3 


xK ( sse DG 45) 


We place the zeros at s = —4 + 2j in order to pull the locus to the left in the s-plane. Comparing the coefficients of 


K3 + K 1 


2 2 
s + 8s +20 =s 
+ + K; K; 


leads to 


and 


— = 20. 
K3 


Therefore Ky = 0.35 and K3 = 0.05. We can now plot a root locus with K K3 as the parameter, as shown in Fig. 10.4. The 
characteristic equation, Eq. (10.5), is 
s? + 8s +20 ) B 


AR (= Do 45) 


The roots for the selected gain, K K3 = 12, lie in the performance region, as shown in Fig. 10.4. The final gains are 


K = 240, 
Kı = |, 
Kə = 0.35, 
K3 = 0.05. 
The controller design results in a settling time of about 1.8 seconds and an overshoot of 3%, as shown in Fig. 10.5. e 


Example 10.3 State-Space System Simulation using LabVIEW 


Consider the system given in state variable form 
X = Ax+ Bu, (10.6) 


epulo = 0) 


In the LabVIEW simulation, you will want to utilize full-state feedback, which means that the entire state vector (in this case, x; and 
X2) needs to be available for feedback. One way to accomplish this is to choose the output matrices C and D such that y — x, thus 


where 


y = Cx + Du. 
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Figure 10.5: Step response for the automatic test system. 
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where 


The full-state feedback control law is 


u = —Kx + v, 


where v is a reference input. The LabVIEW block diagram of the closed-loop system is shown in Fig. 10.6. 


The reference input is supplied by the Sine Wave function. The outputs of the state-space block are the two states x; and x2. 
These states are fed back to the feedback gain matrix block to obtain the feedback signal Kx. The reference input signal v is then 
summed with the (negative) feedback signal forming the plant input signal u = —Kx + v. You can adjust the sine wave signal, 
as shown in Fig. 10.7. The gain matrix K and the system matrices A, B, C, and D are set as illustrated in Figs. 10.8 and 10.9, 
respectively. The feedback gain matrix is chosen as 


K = [1 1]. 
This places the closed-loop poles (that is, the eigenvalues of A — BK at sı = —1 and s2 = —3. Using LabVIEW allows you easy 


access to the feedback gain matrix so that simulation studies using varying feedback gain sets can be readily performed. Notice in 
Fig. 10.9 that the initial conditions are inputs. In the simulation presented here, the initial conditions are chosen to be 


(00) (1 
w= (5 )- (o) 
Once you have defined all the parameters, you can initiate the simulation. The resulting state variable time histories are shown in 


Fig. 10.10. You can use the LabVIEW simulation to conduct numerical experiments easily and to analyze the closed-loop system 
performance under a variety of conditions. o 
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Figure 10.7: Selecting the sine wave parameters on the Sine Wave function. 
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Figure 10.8: Incorporating the feedback gain into the simulation block diagram. 


86 CHAPTER 10 DESIGN OF STATE VARIABLE FEEDBACK SYSTEMS 











Control & Simulation Loop 


M Error + 


Sine Signal Summation State-Space Response Graph 


ee ee - 
— "arcum 
TA State-Space Configuration E + | 


Polymorphic instance Feedthrough Parameter Information 
MIMO i+ | Indirect | + Parameter source — — < 
Parameters s Configuration Dialog Box + | 
Parameter Name Value 
PT SmteSpxe a 
— COS || E initial state (x0) — [10] 
^ El reset? False 
L-Bl reset state (xr) [0] 








» tu Collector 


rnm 








[R][E] 
[Clio] 


patih 


a o 


jasini 









































Figure 10.9: Defining the state variable model via the matrices A, B, C, and D. 
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Figure 10.10: The result of the simulation showing the state variable history. 
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Robust Control Systems 


In this chapter, we investigate robust control systems using LabVIEW. In particular we will consider the commonly used PID 
controller in the feedback control system shown in Fig. 11.1. Notice that the system has a prefilter G p(s). The PID controller has 
the form 
K ps? + Kps+ Kr; 

: i 
The PID controller is not a proper rational function (i.e., the degree of the numerator polynomial is greater than the degree of the 
denominator polynomial). The objective is to choose the parameters Kp, Kp, and Kz to meet the performance specifications and 
have desirable robustness properties. Unfortunately it is not immediately clear how to choose the parameters in the PID controller to 
obtain certain robustness characteristics. An illustrative example will show that it is possible to choose the parameters iteratively and 
verify the robustness by simulation. Using LabVIEW helps in this process, since the entire design and simulation can be mechanized 
utilizing scripts and can easily be executed repeatedly. 


Gc(s) = 


Example 11.1 Robust Control of Temperature 
Consider the feedback control system in Fig. 11.1, where 
G(s) — Gia 
the nominal value is co = 1, and Gp(s) = 1. We will design a compensator based on cg = 1 and check robustness by simulation. 
Our design specifications are as follows: 
1. settling time (2% criterion) T, < 0.5 second, and 


2. minimize the overshoot for a step input. 


For this design we will not utilize a prefilter to meet specification (2) but will instead show that acceptable performance (1.e., low 
overshoot) can be obtained by increasing a cascade gain. The closed-loop transfer function is 


Kps? + Kps + Kr 


LS) ee XUI I NT TER 
s? -F(24 Kp)s-+14+Kp)s+ Kj 


(11.1) 
The associated root locus equation is 
^5? +as +b 


La a V 


Ts) 
Disturbance 
R (s) G(s) Y(s) 
Input li Output 





Figure 11.1: Closed-loop system. 
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Figure 11.2: Root locus for the PID-compensated temperature controller as K varies. 








where 
5 1+K K 
K = Kp +2, E P M 
2+ Kp 2+ Kp 
The settling time requirement T; < 0.5 second leads us to choose the roots of (s? + as + D) to the left of the s = —Cw, = —8 line 


in the s-plane, as shown in Fig. 11.2. To ensure that the locus travels into the required s-plane region, we can choose a — 16 and 
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Figure 11.3: Step response of the PID temperature controller. 


b = 70. We select a point on the root locus in the performance region and we find the associated gain and the associated value of 
c. For our chosen point we find that K = 118. Then with K, a, and b, we can solve for the PID coefficients as follows: 


Kp = K-2=116, 
Kp = a(2+K3)—1=1187, 
Kı = b(2+K3) = 8260. 


To meet the overshoot performance requirements for a step input, we utilize a cascade gain K that will be chosen by iterative methods 
using the CD Step function, as illustrated in Fig. 11.3. The step response corresponding to K = 5 has an acceptable overshoot of 
2%. With the addition of the gain K = 5, the final PID controller is 


Kps? +K K 116s? + 1187s + 8260 
CERK DS” + K ps + | js S op ss | 
S S 


(11.2) 


Now we can consider the question of robustness to changes in the plant parameter co. The investigation into the robustness of the 
design consists of a step response analysis using the PID controller given in Eq. (11.2) for a range of plant parameter variations of 
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Figure 11.4: Robust PID controller analysis with variations in co. 


0.1 < co < 10. The results are displayed in Fig. 11.4. The simulation results indicate that the PID design is robust with respect 
to changes in co. The differences in the step responses for 0.1 < co < 10 are barely discernible on the plot. If the results showed 
otherwise, it would be possible to iterate on the design until an acceptable performance was achieved.The interactive capability of 
LabVIEW allows us to check the robustness by simulation. Q 
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Digital Control Systems 


Many of the LabVIEW functions for continuous-time control design have equivalent counterparts for sampled-data systems. Model 
conversion can be accomplished with the functions CD Convert Continuous to Discrete and CD Convert Discrete to Continuous. 
The function CD Convert Continuous to Discrete converts continuous-time systems to discrete-time systems. Similarly, the 
function CD Convert Discrete to Continuous converts discrete-time systems to continuous-time systems. For example, consider 
the plant transfer function 


Gp(s) = GED 


as shown in Fig. 12.1. For a sampling period of T = 1 second, we determine that 


_ 0.3678(z + 0.7189) 0.3678z + 0.2644 


G(r) = disease 12.1 
©) = GI — 0.3680) ~ 2 — 1.3682 + 0.3680 wen) 


We can use LabVIEW to obtain the G(z), as shown in Fig. 12.2. 


1 


s(s + 1) 





y(t) 






Figure 12.1: A closed-loop sampled data system. 


Example 12.1 Discrete System Unit Step Response 


We can use LabVIEW to compute the response y(kT) of a closed-loop sampled data system using the CD Step Response function, 
shown in Fig. 12.3. Consider again the system shown in Fig. 12.1. The closed-loop transfer function is given by 


0.3678z + 0.2644 


G(z) = ——————À. 
d z? — z + 0.6322 

and the associated closed-loop step response is shown in Fig. 12.3. The zero-order hold is modeled by the transfer function Go(s), 
where 

—sT 


l—-—e 
Go(s) = 


In the next example, we consider the root locus of a sampled-data system. 
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{A Discretization.vi Block Diagram 
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Figure 12.2: Using the CD Convert Continuous to Discrete function to convert G(s) = Go(s)G p(s) to G(z). 


Example 12.2 Root Locus of a Digital Control System 


Consider the system shown in Fig. 12.4 where 


0.3678(z + 0.7189) 
CO) Se ee 
(z — 1)(z — 0.3680) 
and the compensator is 
K (z — 0.3678) 
DR)=— ARA 
z + 0.2400 
The parameter K is a variable yet determined. When 
0.3678(z + 0.7189 
G(z)D(z) = K l ) (12.2) 


(z — 1)(z + 0.2400) ' 


we have the problem in a form for which the root locus method is directly applicable. The rlocus function works for discrete-time 
systems in the same way as for continuous- time systems. Using a LabVIEW script, the root locus associated with Eq. (12.2) is 
easily generated, as shown in Fig. 12.5. 

Remember that the stability region is defined by the unit circle in the complex plane. The LabVIEW function rlocfind can be 
used with the discrete-time system root locus in exactly the same way as for continuous-time systems to determine the value of the 
system gain associated with any point on the locus. Using rlocfind we determine that K = 4.639 places the roots on the unit circle. 
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Figure 12.3: The CD Step Response function generates the output y(kT) for a step input. 
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Figure 12.4: (a) Feedback control system with digital compensation. (b) Block diagram model with G(z) = Z [Go (s)G (s) ]. 
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Figure 12.5: The CD Root Locus function for sampled data systems. 
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